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We extend the previously observed scaling equation connecting the internode distances 
and nodes' degrees onto the case of weighted networks. We show that the scaling takes 
a similar form in the empirical data obtained from networks characterized by different 
relations between node's strength and its degree. In the case of explicit equation for s{k) 
(e.g. linear or scale- free), the new coefficients of scaling equation can be easily obtained. 
We support our analysis with numerical simulations for Erdos-Renyi random graphs with 
different weight distributions. 
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1. Introduction 

The are different ways to represent the structure of real-world networks using the 
notions of graph theory. In particular, depending on the fact if we are interested 
only in the question of existence of connections between nodes or we would also like 
to consider the "intensity" of those links, the described network can be presented 
as unweighted or weighted. In most cases the initial network is of weighted character 
but in simple analytical methods the weights are usually removed. 

The modeling of weighted networks has begun almost parallel with the un- 
weighted case however it is the higher availability of the empirical data that 
made possible to perform extensive and well set research on weighted transporta- 
tion networks, networks of scientific collaboration ^ or in general the networks of 
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social contacts^. I n th e same time the works on connecti ons between topology and 
weight's dynamics ^^El^ transport on weighted networks ^SiZI and optimal paths ^ 
were being developed. Only recently the weighted networks approach has been used 
to create a unified statistics combing Bose-Einstein and Fermi-Dirac distributions 

ini. 

The considerations inclu ded iii this paper are an attempt for extending the 
scaling relations presented in tl^HI] onto the case of weighted networks. In the cited 
works it is shown that the mean distance (or internode distance) between nodes 
with degrees ki and kj is given by the following equation 

(kj) = a - bloghkj, (1) 

where a simple model of unweighted network with a negligible clustering coefficient 
leads to 

log(iV(fc)) , , 1 
a = 1 + 7 ^ " and b=- . (2) 

log K log K 

with being N the size of the network, (fc) its average degree and k the average 
branching degree. The scaling relation ([1} has been spotted in many different net- 
works ranging from public transport networks to biological oneJ^^EI although the 
observed values of the coefficients a and b are different from the results given by 
Eq. ©. 

2. Empirical data 

In the case of weighted network instead of node's degree ki one usually considers 
node's strength Si i.e. 

Si=Y^ Wij (3) 

where Wij is the weight of the link between nodes i and j. In order to find the 
relation between the internode distances {Itj) and the product of nodes' strengths we 
examined the following systems: the network of flights between world-wide airport^, 
public transport network in Warsawj and two networks made from connections 
between the web-pages of web portalcl- In the case of the airport network the weight 
Wij stands for the number of different carriers making the connections between 
the airports i and j while the weights in Warsaw public transport network are 

''The data contained information about N = 3073 airports with average degree value (k) = 10.3 
and they were obtained from the web- page of OAG http://www.oag.com corresponding to the 
status on 14th March 2006. 

''The network is described in 1121 with A'^ = 1530 stops and average degree N = 1530 (fc) = 2.8; 
data was downloaded form Warsaw Tran spor tation Authorities web-page in April 2004. 
"^The data from two large Polish portals were obtained thanks to collaboration with Gemius 
company. For reasons of information confidentiality they will be further referred to as portal A 
(number of nodes N = 195, average degree (fc) = 70) and portal B (N = 512, (fc) = 102) 
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proportional to the total number of buses and trams traveling between stops i and 
j. In the web portals data the weight was defined as the total number of users 
who, being first at page i surfed to page j. For simplicity reasons in all datasets 
the weights were symmetrized by calculating the average of the original values i.e. 
w = {wij + Wji)/2. Due to high reciprocity between Wij and Wji this does not affect 
further analysis. 




Fig. 1. The average node's strength versus node's degree {s{k)) for world-wide airport network 
(A) and pubhc transport network in Warsaw (C). The solid line is a fit to the empirical data 
{s{k)) ~ (A) and {s{k)) = Sk — 6 (C). The dependence of the internode distances {lij) on 

the product of nodes' strengths SiSj (circles) and on the product of nodes' degrees kikj rescaled 
by the relation {s{k)) (triangles) for the airport connections (B) and public transport in Warsaw 
(D). The solid line in (B) is defined by the rescaled coefficients a and b obtained from Eqs II14I I. 

The analysis of the results obtained from the above described datasets is pre- 
sented in Figures [T] and [5] For each network we plot the dependence of average 
node's strength and its degree (s(fc)) (Figures [TJ\, [Tp, [DA. and[2P) as well as the 
internode distances {hj) separately as a function of the product of degrees kikj and 
as a function of the product of strengths SiSj (Figures [T}3,[Tp, 03 and[2P). 

The plots presenting (/y) as a function of kikj were rescaled with respect to 
node degree values, substituting it with an approximated relation between the node 
strength and its degree, marked as solid curve in {s{k)) plots. The {s{k)) equa- 
tion can take different forms: starting from lack of correlations in case of public 
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transport network in Warsaw ((s(fc) ^ k), through a scale- free relation in airport 
network ((s(fc)) ~ k"), ending at a complicated dependence seen in web portals 
data ((s(fc)) ~ k" exp{/3k)). Nevertheless in each example the coincidence between 
the rescaled data and the empirically observed relation {lij{siSj))is clearly visible. 
The oscillations observed on the logarithmic trend seen in Figures [2j3 and [2jl) are 
an effect of the high average degree present in both portals data. The origin of this 
phenomenon is extensively described in^i^l 




Fig. 2. The average node's strength versus node's degree {s{k)) for Portal A (A) and portal 
B (C). The solid line is a fit to the empirical data (s(fc)) = 2fc^'^ exp(fc/30) (A) and {s(fc)) = 
0.5fc^ exp(fc/100) (C). The dependence of the internode distances (lij) on the product of nodes' 
strengths SiSj (circles) and on the product of nodes' degrees kikj rescaled by the relation (s(fc)) 
(triangles) for the Portal A (B) and Portal B (D). 



3. Scaling in the absence of s(fc) correlations 

The scaling relation observed in the absence of strength-degree correlations can 
be easily explained using the approach previously applied for internode distances 
jj^Iini fact, considering a branching process one obtains the following equation 
connecting the degrees of nodes i and j, the distance x between them, the size of 
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the network TV, its average degree (fc) and the average branching degree k: 

= ^. (4) 

Multiplying the nominator and denominator at both sizes of the equation by the 
average value of the weight in the network {w) and assuming that in the absence of 
the s{k) correlations the following relation is satisfied 

Si = (w)fci, (5) 

one gets 

N{k){wf 

SiK = . (6) 



This equation immediately results in the form of the internode distances scaling 
relation as a function of nodes' strengths 

{kj) ^a-b\ogSiSj, (7) 

where coefficients a and b are given by 

^ log(iV(fc)(w)2) ~ 1 

log K log K 

The above considerations lead to the following conclusions: (a) the slope b of the 
scaling function is the same as in the unweighted case, (b) the relation between 
coefficients a and a takes the form 

log(w)^ 

a = a-\ . (9) 

logK 

In order to compare the consistency of the Eqs ([7]) and ^ with the numerical 
simulations we performed the trails (see Figure [3|) for the Erdos-Renyi random 
network (ER network) of iV = 1000 nodes and (fc) — 5. The weights were randomly 
drawn from the Gauss distribution 

p^H = -=^e (10) 
V 27r(T,„ 



characterized by the mean value (w) = 10 and three different values of the standard 
deviation: cr^ = 0.1 (Figures[5J^ and[2P), cr„ = 1 (Figures[3j3 and[3|5) and a^, = 3 
(Figures [3p and[5p). A comparison of those weight distributions is presented in 
Figure [3fl. As one can see the width of the distribution plays the key role with 
respect to the shape of the strength distribution P{s); in case when (t„, = one 
gets of course the Poisson distribution 

P(.)^e-W-i^, (11) 
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for s — (w) ,2{w) ,3{w) ^ .... In general, the strength distribution takes the form 

oo / k \ 

P{s) = J2 P(fc)Prob E = ' (12) 

fe=l \i=l / 

assuming that weights Wi are i.i.d. random variables drawn from the distribution 
Pwiw). Taking into account the considered distributions p{k) and Pw{w) one arrives 
at the relation 

P(,) ^ y ^^e-^^e-<^> (13) 

The strength distributions taken from the numerical simulations for three different 
values of as well as the values P{s) numerically summed using Eq. (jl3p are shown 
in Figures [3j^, |3j3 and[3lC The characteristic oscillations seen in P{s) distribution 
around the multiply of the average weight ({w) ,2{w) ,3{w) ,...) vanish along with 
the increasing value of which in result gives the distribution a smooth shape. 
The form of the strength distribution changes the relation between the internode 
distances and the product of strengths (see Figures [3jI)-F). The relation bewteen 
the average node's strength and node's degree follows Eq. ([5]) and is shown in Figure 

For (Tu, = 0.1 (Figure [3jl)) the points obtained from the numerical simulations 
overlap those predicted by Eq. ^ with the coefficients a and b defined by the 
parameters of the ER network and the average weight (w) . A comparison between 
relations {lij){siSj) and {lij)(kikj) is presented in Figure |3j showing that indeed 
the slopes in both cases are identical as predicted by Eqs ([2]) and ([8|. The increase 
in the value leads to emergence of the oscillations on the log{siSj) trend line 
(Figure [Sfl) and eventually for large values of the scaling relations breaks down. 
This effect is illustrated in Figure |3p. 



4. Scaling in the presence of s(fc) correlation 

Section [2] contains examples of the scaling of internode distances on the product of 
nodes' strengths in the presence of non-linear s(k) correlations. Previously proposed 
branching-tree model has no application in this case, as one should sum up all the 
weights on the tree's surface, which is impossible for any general form of s{k). 
Certainly, if s(fc) has an explicit form, there is a possibility of reversing the relation, 
putting it into Eq. ([4]) and obtaining the coefficients 5 and b defining the {Uj) relation 
as a function of siSj. For instance, if s{k) = ^fc", reversing gives k = (s/A)^/" and 
as a result the rescaled coefficients a and b will take the form 

6 = - and a = a + 2 log Ab. (14) 
a 

Such an example is presented in Figure [T)3 (solid line) where, after getting coeffi- 
cients a and b from fitting to relation ^ they were put into Eqs I14p thus obtaining 
the values a and b. 
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Fig. 3. Numerical simulations of the ER network with N = 1000 and (fc) = 5 and the weights 
coming from Gauss distribution llOI I with (to) = 10 and different values of standard deviation a^. 
(A-C) Strength distributions for = 0.1 (A), CTu, = 1 (B) and = 3 (C); solid lines correspond 
to numerically summed Eq. I I13I I. (D-F) Internode distances versus the product of strengths for 
a-w =0.1 (D), (Jm = 1 (E) and o-^, = 3 (F); solid lines correspond to Eq. JT]!. (G) The node's 
strength versus its degree; solid line corresponds to equation s{k) = 10k. (H) Weight distributions 
for (Tu, = 0.1 (circles), = 1 (triangles) and cr^ = 3 (squares); solid lines come from Eq. mOI I. (I) 
A comparison between relations {lij){siSj) (circles) and {lij){kikj) (triangles) for = 0.1; solid 
lines come from Eqs ((1} and 10. 

For more systematical study we performed numerical simulations for ER net- 
work with N = 1000 nodes and {k) — 5, with specific s(k) relations. To generate 
the correlations between node's strength and its degree we used the observations 
from work : if in the given network there are no degree-degree correlations (the 
asortativity coefficient r is equal to 0) and the edge between nodes i and j has 
weight (wij), proportional to the expression (kikj)" then the strength of node i is 

~ h{w,,) ^ fci+". (15) 

Using the above described procedure, we imposed weights onto ER network accord- 
ing to relations s{k) = ISOfc^ (Figure H^) and s{k) = 80A: sin(7rfc/18) (Figure gp). 
As one can observe in Figures [4j3 and|4jD the expected overlap of relation {kj) as 
a function of SiSj (circles) with the rescaled data of {lij)(kikj) (triangles) was ob- 
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tained. Moreover, in Figure |4j3 there can be seen a good fitting of data with linear 
function determined by the coefhcients taken from the Eqs (|14p. 




Fig. 4. Plots A and C: the average node's strength versus its degree k. Plots B and D: the 
internode distance versus the product of nodes' degrees kikj rescaled by the relation between 
node's strength and its degree (triangles); the internode distance versus the product of nodes' 
strengths SiSj (circles). Plots A and B were obtained for ER network with N = 1000, (fc) = 5 
where weights fulfill the equation Wij ~ (kikj)^ , while in the case of plots C and D a corresponding 
relation was Wij ~ sin(fci) sin(fcj). 



5. Conclusions 

The relation between internode distances and the product of the nodes' strengths 
is an extension of the previously examined dependence linking (lij) and kikj. As 
expected the scaling takes a similar form while the coefficients a and b are altered by 
the relation between average node's strength and its degree. In the case of explicit 
equation for s{k) (e.g. linear or scale-free), the new coefficients a and b can be easily 
obtained. A breakdown in scaling relation while increasing the width of the weight 
distribution has also been observed. 
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